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1 Introduction 



In this paper, we study a discrete-time approximation for the solution of a forward-backward 
stochastic differential equation (FBSDE) with a jump and taking the following form 



X t = x+ f b{s,X s )ds+ [ a(s)dW s + [ (3{s,X s -)dH s , 
Jo Jo Jo 

Y t = g(X T )+[ f( Sj X Sj Y Sj Z s ,U s )ds- [ Z s dW s - [ U s dH s , 



< t < T, 



where H t = l T <t an d t is a jump time, depending on W, which can represent a default 
time in credit risk or counterparty risk. Such equations naturally appear in finance, see 
for example Bielecki and Jeanblanc [2], Lim and Quenez [19], Ankirchner et al. |T] for an 
application to the exponential utility maximization problem and Kharroubi and Lim [15] for 
the hedging problem in a complete market. The approximation of such equation is therefore 
of important interest for practical applications in finance. 

This paper is the second part of a two sides work studying the discretization of such 
FBSDEs. The first part [16] deals with the case of a Lipschitz continuous generator / and 
this second part focuses on the case of a generator / with quadratic growth w.r.t. z. 

In the literature, the problem of discretization of FBSDEs with Lipschitz generator has 
been widely studied in the Brownian framework, i.e. no jump, see e.g. [201 d El El EH |6]. 
More recently, the case of quadratic generators w.r.t. z has been considered by Imkeller et 
al. [8] and Richou [23] . 

The discretization of BSDEs with jumps has been studied by Bouchard and Elie [I] in 
the case where the generator is Lipschitz continuous and the jumps are independent of the 
Brownian motion. 

In this paper, we study the discrete-time approximation of FBSDEs with 

- a generator satisfying a quadratic growth w.r.t. the variable z, 

- a jump time r that could be dependent on the Brownian motion via the density 
assumption. 

Such FBSDEs are of interest in finance since utility maximization problems lead in 
general to quadratic generators and default time models generally impose the dependence 
of r w.r.t. the Brownian motion. 

To study the discretization of such equations, we can not directly work on the BSDE as 
done in [I] since their approach uses a regularity result for the process Z based on Malliavin 
calculus. In our situation, the dependence of r w.r.t. W prevent us from using such an 
approach since no Malliavin calculus theory has been set for this framework. 

To get a discrete-time approximation scheme we then use the results of [T5], that allow 
to decompose the FBSDE with a jump in a recursive system of Brownian FBSDEs. We 
then provide estimates on the solutions to each Brownian FBSDE. These estimates allow 
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to prove that the solution satisfies a Lipschitz FBSDE. Finally, we obtain a discrete-time 
approximation scheme by using the results of the previous part [16] on Lipschitz BSDEs. 

The paper is organized as follows. The next section presents the FBSDE, the different 
assumptions on the coefficients of the FBSDE and recalls the result of [T5]. In Section 3, we 
give some estimations on the solution of the FBSDE. In Section 4, we give a discrete-time 
approximation scheme for the FBSDE and provide a global error estimate. 

2 Preliminaries 
2.1 Notation 

Throughout this paper, we let (Q,Q,¥) a complete probability space on which is defined a 
standard one dimensional Brownian motion W. We denote F := (^t)t>o the natural filtration 
of W augmented by all the P-null sets. We also consider on this space a random time r, 
i.e. a nonnegative J-"-measurable random variable, and we denote classically the associated 
jump process by H which is given by 

H t := t T <t, t>0. 

We denote by D := (T> t ) t > the smallest right-continuous filtration for which r is a stopping 
time. The global information is then defined by the progressive enlargement G := (Qt) t >o °f 
the initial filtration where 

Gt := p| (j- t+£ v% £ ) 

for all t > 0. This kind of enlargement was introduced by Jacod, Jeulin and Yor in the 80s 
(see e.g. [TO], [IT] and [S]). We introduce some notations used throughout the paper 

- V(F) (resp. V(G)) is the cr-algebra of F (resp. G)-predictable measurable subsets 
of Q x M + , i.e. the a-algebra generated by the left-continuous F (resp. G)-adapted 
processes, 

- VA4(F) (resp. VA4(G)) is the a-algebra of F (resp. G)-progressively measurable 
subsets of Q x M + . 

We shall make, throughout the sequel, the standing assumption in the progressive enlarge- 
ment of nitrations known as density assumption (see e.g. [121 H2J HSj ) - 

(DH) There exists a positive and bounded V(¥) eg) £>(lR + )-measurable process 7 such that 

P[r G d6 I T t ] = -y t (6)d6 , t > . 

Using Proposition 2.1 in [15] we get that (DH) ensures that the process H admits an 
intensity. 
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Proposition 2.1. The process H admits a compensator of the form Xtdt, where the process 
X is defined by 



-t t<T , t > . 



P[r >t | J=i] 

We impose the following assumption to the process A. 
(HBI) The process A is bounded. 

We also introduce the martingale invariance assumption known as the (H)-hypothesis. 
(H) Any F-martingale remains a G-martingale. 

We now introduce the following spaces, where a,b G R+ with a < b, and T < oo is the 
terminal time. 

- Sfe[a,b] (resp. <Sf°[a, &]) is the set of VAi(G) (resp. T^M (F))-measurable processes 
(Xt)te[a,b] essentially bounded 

ll*1ls°°ra&i := esssu Pl^*l < 00 • 

te[a,b] 

- <Sg[a,6] (resp. S^[a,b]), with p > 2, is the set of VA4(G) (resp. VM. (F))-measurable 
processes (5^)te[ ,&] suc h that 



\Y\ 



SP[a,b] 



E 



sup \Y t \ 

L t£[a,b] 



< OO 



Hq[ci, b] (resp. ii/|[a, 6]), with p > 2, is the set of 7 7 (G) (resp. "P(F))-measurable 
processes (Z t )te[a,b] sucn that 

r>6 



l#P[a,&] 



E 



< 00 



- £ 2 (A) is the set of P(G)-measurable processes (U t )t£[o,T] such that 



\U\ 



E 



XJUJ ds 



'-jo 



< 00 . 



2.2 Quadratic growth Forward-Backward SDE with a jump 

Given measurable functions b : [0, T] x R -> E, a : [0, T] -> R, /3 : [0, T] x R -> E, 5 : K -> R 
and / : [0, T] x E x R x R x R — > R, and an initial condition x G R, we study the discrete-time 
approximation of the solution {X, Y, Z, U) in «S£[0, T] x <Sg°[0, T] x H&[0, T] x L 2 (A) to the 
following forward-backward stochastic differential equation 

X t = x+ f b{s,X s )ds+ [ a(s)dW s + [ /3(s,X s -)dH s , 0<t<T, (2.1) 
Jo Jo Jo 

Y t = 9(X T ) + £ f(s,X s ,Y s ,Z s ,U s (l-H s ))ds 

(2.2) 



Z s dW s 



U s dH s , < t < T 



when the generator of the BSDE (12. 2p has a quadratic growth w.r.t. z. 



Remark 2.1. In the BSDE ( 12. 2ft . the jump component U of the unknown (Y, Z, U) appears 
in the generator / with an additional multiplicative term 1 — H. This ensures the equation 
to be well posed in «Sg°[0, T] x H^[0,T] x L 2 (X). Indeed, the component U lives in L 2 (X), 
thus its value on (r A T,T] is not defined, since the intensity A vanishes on (r A T, T]. We 
therefore introduce the term 1 — H to kill the value of U on (r A T, T] and hence to avoid 
making the equation depends on it. 

We first prove that the decoupled system (I2.ip - (I2.2I) admits a solution. To this end, we 
introduce several assumptions on the coefficients b, a, (3, g and /. We consider the following 
assumptions for the forward coefficients. 

(HF) There exist two constants K a and L a such that the functions b, a and /3 satisfy 

\b(t,0)\ + \a(t)\ + \l3(t,0)\ < K a , 

and 

\b(t,x)-b(t,x')\ + \P(t,x)- P(t,x')\ < L a \x-x'\, 
for all (t, x, x') G [0,T] xlx R. 

For the backward coefficients g and /, we consider the following assumptions. 
(HBQ) 

- There exist three constants M g , K g and K q such that the functions g and / satisfy 

l<?0)l < M 91 

\g(x) - g(x')\ < K g \x - x'\ , 

\f(t,x,y,z,u)-f(t,y',z,u)\ < K q \y-y'\, 

\f(t,x,y,z,u)\ < K q {\ + \y\ + \z\ 2 + \u\) , 

for all (t, x, x\ y, y\ z, u) G [0, T] x R 2 x R 2 x R x R. 

- For any R > there exists a function mc^ R such that lim^o mc^(e) = and 

\f(t,x,y,z,u-y)- f(t,x,y',z',u-y')\ < mc f R (e) 

for all (t,x,y,y',z,z f ,u) E [0, T] x R x R 2 x R 2 x R s.t. \y\, \z\, \y'\, \z'\ < R and 
\y — y'\ + \z — z'\ < e. 

- f(t,.,u) = f(t,.,0) for all m G R and all t G (r AT,T}. 

- The function f(t,x,y,.,u) is convexe (or concave) uniformly in (t,x,y,u) G [0,T] x 
R x R x R. 
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In the sequel K denotes a generic constant appearing in (HBQ) and (HF) and which may 
vary from line to line. 



In the purpose to prove the existence of a solution to the FBSDE (|2.1[) - (j2.2p we follow the 
decomposition approach initiated by [T3] and for that we introduce the recursive system of 
FBSDEs associated with f[2~ll- f[2~2p . 

• Find (X 1 (9),Y 1 (e),Z 1 (e)) G «S|[0,T] x «S£°[0,T] x H$[6,T] such that 

Xft9) = x+ f b(s,X s L (9))ds + [ a{s)dW s + (5{e,Xl-(6))t e <t , 0<t<T, (2.3) 
Jo Jo 

Y t \9) = g{X 1 T {9))+J^ f(s,X 1 s (6)X(0),Z 1 s (e),0)ds- J Z]{9)dW s , 9 < t < T , (2.4) 
for all 9 G [0,T]. 

• Find (X°,Y°,Z ) G <Sf[0,T] x <Sf°[0,T] x H§[0,T] such that 

X° = x+ [ b( Sl X° s )ds+ [ a(s)dW s , 0<t<T, (2.5) 
Jo Jo 

y? = 9{X°t)+ [ f(a,X° a ,Y?,Z°.X(s)-Y?)ds- jf Z° a dW a , < t < T . (2.6) 

Then, the link between the FBSDE (l2~T]) -(l2T2l and the recursive system of FBSDEs Q - 
(12 .4p and (I2.5p - (I2.6I) is given by the following result. 

Theorem 2.1. Assume that (DH), (HBI), (H), (HF) and (HBQ) hold true. Then, the 
FBSDE fPO)- (l2~2p admifo a unigue sotozon (X, Y, Z, U) G S£[0, T] x 5^[0, T] x H^[0, T] x 
L 2 (A) gwen by 

{ X t = XX<r + Xl(r)l T < t , 
Y t = Y?t t<T + Y?{T)t T < t , 
Z t = ZX<T + Zl(r)l T<t , 
( U t = (Y t \t)-Y t °)t t < T , 

where (X 1 (9), Y 1 (9), Z 1 (6)) is the unique solution to the FBSDE Q-O in «S|[0,T] x 
S£°[9,T] x H$[8,T], for 9 G [0,T], and (X°,F°,Z ) is the unique solution to the FBSDE 
([23]) -fl2D in SI [0, T] x «S|° [0, T] x fl|[0, T] . 

Proof. The existence and uniqueness of the forward process X and its link with X° and 
X 1 have already been proved in the first part of this work [16]. We now concentrate on the 
backward equation. 

To follow the decomposition approach initiated by the authors in [T3] , we need the generator 
to be predictable. To this end, we notice that in the BSDE 02.2p . we can replace the generator 
(t,y,z,u) H> f(t,X t ,y,z,u(l-H t )) by the predictable map (t,y,z,u) n- f(t,X t -,y,z,u(l- 
H t -)). 
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Suppose that (DH), (HBI), (H) and (HBQ) hold true. The existence of a solution 
(Y, Z, U) e Sq[0, T] x Hq[0, T] x L 2 (X) is then a direct consequence of Proposition 3.1 in [15] . 
We then notice that from the definition of H we have f(t, x, y, z, u(l — H t )) = f(t, x, y, z, 0) 
for all t e (r A T,T\. This property and (DH), (HBI), (H) and (HBQ) allow to apply 
Theorem 4.2 in [T5], which gives the uniqueness of a solution to BSDE (12. 2p . □ 

Throughout the sequel, we give an approximation of the solution to the FBSDE ( 12. ip - 
( 12. 2 p by studying the approximation of the solutions to the recursive system of FBSDEs 
Q-flHD and ( I2~5|) - (l23|) . 

3 A priori estimation on the gain process 

Before giving the discrete-time scheme for the FBSDE ( I2.ip - ()2.2p we give a uniform bound 
for the processes Z° and Z 1 which allows to prove that the BSDE ( I2.2p is Lipschitz and 
thus we can use the discrete-time scheme given in [16]. For that we introduce the BMO- 
martingales class, and we also give some bounds for the processes X°, X 1 , Y° and Y 1 . 

3.1 BMO property for the solution of the BSDE 

To obtain a uniform bound for the processes Z° and Z 1 we need the following assumption. 
(HBQD) There exists a constant Kf such that the function / satisfies 

\f(t : x : y,z : u)-f(t',x' : y',z',u')\ < K f [\x-x'\ + \y-y'\ + \u-u'\ + \t-t'\ 1 *} 

+£/,*(! + M + \z'\)\z-z'\ , 

for all (t, x, x', y, y\ z, z' : it, u') G [0, T] x I 2 x I 2 x I 2 x M 2 . 



In the sequel of this section, the space of BMO martingales plays a key role for the a 
priori estimates of processes Z° and Z l . We refer to [2] for the theory of BMO martingales. 
Here, we just give the definition of a BMO martingale and recall a property that we use in 
the sequel. 

Definition 3.1. A process M is said to be a BMOf[0,T] -martingale if M is a square 
integrable ¥ -martingale s.t. 



\M\ 



BMO V [0,T] 



sup E 

reT F [0,T] 



M T - Mr 



1/2 



< oo 



where 7f[0, T] denotes the set of ¥ -stopping times valued in [0, T]. 

The BMO condition provides a property on the Dolean-Dade exponential of the process M. 



S 



Lemma 3.1. Let M be a BMOf[0,T]-martingale. Then the stochastic exponential S(M) 
defined by 

£(M) t = exp(M t -i<M,M)) 0<t<T, 

is a uniformly integrable ¥ -martingale. 
We refer to [H] for the proof of this result. 



We first state a BMO property for the processes Z° and Z 1 , which will be used in the 
sequel to provide an estimate for these processes. 

Lemma 3.2. Under (HF), (HBQ) and (HBQD), the martingales J Z Q s dW s and J Z](9)t s > e dW s , 
9 G [0, T] are BMOf[0,T]-martingales and there exists a constant K which is independent 
from 9 such that 



Z" s dW s 



sup 

0e[o,T] 

Proof. Define the function q 

(j)(x) 

We notice that 6 satisfies 



Zl(9)t s > e dW s 



->• R by 



BMO F [0,T] 



BMO F [0,T] 



< K 



< K. 



r £ 2K q x 



-2K g x-l)/2K* xe 



(3.1) 



(j)'(x) > and -<f>"(x) - K q <j/{x) = 1 , 

for x > 0. Since Y° and Y 1 ^) are solutions to quadratic BSDEs with bounded terminal 
conditions, we get from Proposition 2.1 in [18] the existence of a constant m such that 

(3.2) 



|^ ||s°°[o,t] < rn and sup ||y 1 (6')|| S o O [ 0)T ] < m 

6e[o,T\ 



Applying Ito's formula we get 

r-T 



0(i;° + m)+E(^ ^"(Y s ° + m)\Z s fds\^) = 

E(0(F T ° + m) \JF V ) + <t>'(Y s ° + m)f(s, X° s , Y s °, Z° s , Y s \s) - Y s °)ds 

for any F-stopping time v valued in [0,T]. From the growth assumption on the generator / 
in (HBD), (ED and 1ET2]) . we obtain 

<t>{Y° + m)+E(^ \Z° s \ 2 ds\^ < 

E(0(F T ° + m)|^)+E( / ^(Y s ° + m)K q (l + 2\\Y°\\ soc + sup \\Y 1 (9)\\ s ^ [&> T\)ds T u 

v Jv ee\o,T] 
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This last inequality and (13. 2p imply that there exists a constant K which depends only on 
m, T and K q such that for all F-stopping times v G [0, T] 



E 



T 



\Z°\ 2 ds 



J 7 J < K. 



For the process Z 1 , we use the same technics. Let us fix 9 G [0, T]. Applying Ito's fomula 
we get 



cf>(Y^ e {9)+m)+E 



1 



<P"(Y s \e) + m)\Zm\ 2 ds 



E(4>(Y}(0) + m)\F v )+E( / cf>'(Y s 1 (9) + m)f(s,X 1 s (9),Y s 1 (d),Z 1 s (9),0)ds 



u\/6 



for any F-stopping time v valued in [0,T]. From the growth assumption on the generator / 
in (HBQ), (ED) and (E2D, we obtain 



<P(YXe(0) + m)+E 









( f m 


)\ 2 ds 




v Juve 







T 



+ E( / ( j ) '(Y s 1 (e)+m)K q (l + \\Y 1 (6)\\ Sac[e , T] )ds 



This last inequality and (13. 2\\ imply that there exists a constant K which depends only on 
m, T and K q , such that for all F-stopping times v valued in [0,T] 



\zl(e)\ 2 t s > e ds 



^ K 



□ 



3.2 Some bounds about X° and X 1 

In this part, we give some bounds about the processes X° and X 1 which are used to get a 
uniform bound for the processes Z° and Z . 



Proposition 3.1. Suppose that (HF) holds. Then, we have 



|VX°| := 
and for any 9 G [0,T] we have 



dX» 



dx 



< e LaT , < t < T 



\vxt 



V e Xl(9)\ ■■- 

dX}{9) 



dX}{9) 



< e LaT , 9 <t<T 



dx 



dX}{9) 

< (1 + L a e LaT )e LaT , 9 <t<T 



(3.3) 



(3.4) 



(3.5) 
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Proof. We first suppose that b and f3 are C\ w.r.t. x. By definition we have 



ft 

VX t u = 1+1 d x b{s,X°)VX° s ds , 0<t<T. 



o 



We get from Gronwall's lemma 

|VX t °| < e LaT , < t < T . 

In the same way, we have 

V e Xl(6) = l+[ d x b(s,X 1 s (9))V d X 1 s (9)ds , 9<t<T, 
Je 

and from Gronwall's lemma we get 

IV^ 1 ! < e LaT , 9<t<T. 
Finally we prove the last inequality. By definition 

vx, 1 ^) = i + / d x b( Sj x 1 s (e))vx 1 s (e)ds + d x p{9,x° )vx° e , e<t<r . 

Jo 

Using the inequality (13. 3p . we get 

\VX}(6)\ < l + L a e LaT + [ L a \VXl(6)\ds, 9<t<T, 

Jo 

from Gronwall's lemma we get 

\VX}(9)\ < (1 + L a e LaT )e LaT , 6<t<T . 

When b and (3 are not differentiable, we can also prove the result by regularization. We con- 
sider a density q which is C£° on R with a compact support, and we define an approximation 
of (6,/3) inC^by 

(b e ,(3 e )(t,x) = -f(b,/3)(t,x')q(?—^-)dx', (i, x) e[0,T]xl, 

We then use the convergence of (X 0,e , X 1,e (#)) to (X°,X 1 (6 1 )) and we get the result. □ 

3.3 Some bounds about Y° and Y 1 

In this part, we give some bounds about the processes Y° and Y 1 which are used to get a 
uniform bound for the processes Z° and Z 1 . 

Lemma 3.3. Suppose that (HF), (HBQ) and (HBQD) hold. Then, for any 9 e [0,T] 

dYt ^ d) < e^ +K ^ T (K g + TK f ) , 9<t<T. (3.6) 



dX}(9) 
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Proof. We first suppose that b, f and g are C\ w.r.t. x, y and z. In this case (X 1 (9), Y l (9), Z l {9)) 
is also differentiable w.r.t. Xg(9) and we have 



V e Z l s (9)dW s 

for t G [9, T}. Define the process R{9) by 

Rt{9) := exp^dyf^XMXW^M^s^ds) , 0<t<T. 

Applying Ito's formula, we get 

Rt{9)V e Y t \9) = R T {9)Vg(X l T (9))V e X l T (9) 

+ J R s (9)d x f(s, Xl(9), Y s \9), Z l s (9), 0)V e Xl(9)ds 

- J R s {9)V e Z]{9)dWl{9) , 9<t<T , 

where the process W 1 (9) is defined by 

W t \9):=W t - [ d g f(s 7 X 1 a (e)X(0),ZM,O)*s>eds 



(3.7) 



{3.1 



(3.9) 



for t e [0, T]. From (HBQD), there exists a constant K > such that we have 



d z f(s, Xl(9), Y s 1 (9), Zl(9), 0)l s > e dW s 



BMO F [0,T] 



K[l+ sup E 

i?S7i[0,T] L Jfi 



Z x A9)\H s>e ds 



< 



< 



K 1 



Zl(6)l a > e dW s 



BMO V [0,T] 



< oo 



where the last inequality comes from Lemma 13. 21 

Hence by Lemma O the process S(J Q d z f(s, X}(9), Y^{9), Z\(0), 0)t s > e dW s ) is a uni- 
formly integrable martingale. Therefore, under the probability measure Q 1 ^) defined by 

dQ\9) 



dF 





- °u 


Ti 





e^J d z f(s,Xl(0),Y}{6),Zl{6),O)l a > e dW.) t , < t < T 



we can apply Girsanov's theorem and W 1 (9) is a Brownian motion under the probability 
measure Q 1 ^). We then get from ( 13. 8ft 



Rt{0)V e Y? = E Q1(9) R T (9)Vg(X±)V 6 X± + / R s (9)dJ{s, X\, Y s \ Z], 0)V e X]ds 
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This last equality, (HBQD) and H3.4[) give 

\V e Y^{9)\ < e {La+K f )T (K g + TKf) , 6 <t<T . (3.10) 
When b, f and g are not C^, we can also prove the result by regularization as for Proposition 

o □ 

Lemma 3.4. Suppose that (HF), (HBQ) and (HBQD) hold. Then, 

|Vy/(t)| < (1 + L a e LaT )e (La+K f )T (K g + TKf) , 0<t<T . (3.11) 
Proof. Firstly, we suppose that b, (3, g and / are C\ w.r.t. x, y and z. Then, for any 

te [o,T] 

8Y 1 (t) 

VY'it) := = Vg(X&t))VXUt) (3.12) 

+ jf Vf(s, Xl(t), Y s \t), Z l s (t), 0)(VX s 1 (t), VY'it), VZ]{t))ds 

T 

VZ]{t)dW s . 

Applying Ito's formula, we get 

R t (t)VY t \t) = R T (t)Vg(X^(t))VX^(t) 

+ J R s (t)d x f(s, Xl(t), Y s l (t), Z l s (t), 0)VXl(t)ds 

- J R s {t)VZl{t)dWl{t) , < t < T , (3.13) 

where the process W /1 (.) is defined in (13.91) . We have proved previously that W l {t) is a 
Brownian motion under the probability measure Q : (t). We then get 

Rt{t)VY t \t) = E mt) \R T (t)Vg(X 1 T (t))VX 1 T (t)+ [ R s (t)d x f(s, Xl(t), Y^t), ^(t), 0)VXl{t)ds 



This last inequality, (HBQD) and flUD give 

\VYt(t)\ < (1 + L a e LaT )e {La+K f )T (K g + TK f ) , < t < T . 
When b, f and g are not C^, we can also prove the result by regularization as for Proposition 

o □ 

Lemma 3.5. Suppose that (HF), (HBQ) and (HBQD) hold. Then, 

|Vy t °| < e {2K f +La)T (K g + K f T)(l + TK f e K f T {l + L a e LaT )) , < t < T . 
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Proof. We first suppose that b, (3, g and / are Cl w.r.t. x, y, z and it, then (X°, Y°, Z°) is 
different iable w.r.t. x and we have 



rT 



(3.14) 



+ ^ (vf(s,X s ,Y s ,Z s ,Y s \s) -Y s °)(VX s ,VY s ,VZ° s ,VY s l (s) - VY s °)ds 

VZ° s dW s , < t < T . 
Define the process R° by 



R 



o 



exp ( j\d y - d u )f(s, X° s , Y s °, Z° s , Y s \s) - Y s °)ds) , < t < T . 



Applying Ito's fomula we have 



RtVY t " = R^Vg(X^)VX} 

T 



+ J R»dJ(s, X^ F s u , Zl Y s \s) - Y»)VX»ds 
+ jf R° s d u f(s, X° s , Y s °, Z° s , Y s \s) - Y s °)VY s l (s)ds 



l° s VZ° s dW° 



where dW° := dW s - d z f{s,X° s ,Y^Z° s ,Y}{s) - Y s °)ds. From (HBQD), there exists a 
constant K > such that we have 

d z f(s } X° s , Y s °, Zl Y s \s) - Y s °)dW s "' 



< 



BMOw[0,T] 









(l + sup E 


' [\Zl\ 2 ds 




^ i?S7i[0,T] 







K 1 



Z" s dW s 



BMO V [0,T] 



< 



< oo 



where the last inequality comes from Lemma 13.21 

Hence by Lemma [3TT1 the process £(f d z f(s, X®, Y®, Z®, Y}(s) — Y®)dW s ) is a uniformly 
integrable martingale. Therefore, under the probability measure Q° defined by 

S( ( d z f(s,X s ,Y s ,Z° s ,Y s \s) ~Y s °)dW s 



dF 



ft 



we can apply Girsanov's theorem and W° is a Brownian motion under the probability 
measure Q°. Then, we get 

R° t VY t ° = E Q o Lr° Vg(X%)VX% 

+ £ R%f(s, X°, Y s °, Zl Y s \s) - Y s °)VX° s ds 



+ / R%f(s,X s X,Z° s ,Y s 1 (s)-Y s )VY s 1 (s)ds 
11 



Using inequalities (I3.3P and (13. lip we get 

\VY?\ < e {2K f +La)T (K g + K f T) (l + K f Te K f T (l + L a e LaT )) , < t < T . 

When b, (3, f and g are not C£, we can also prove the result by regularization as for 
Proposition 13.11 □ 

3.4 A uniform bound for Z° and Z 1 

Using the previous bounds, we obtain a uniform bound for the processes Z° and Z l . 

Proposition 3.2. Suppose that (HF), (HBQ) and (HBQD) hold. Then, for any 9 e 
[0, T], there exists a version of Z 1 (9) such that 

\Z\{d)\ < e {2La+K f )T (K g + TK f )K a , 6<t<T . 

Proof. Using Malliavin calculus, we have the classical representation of the process Z l {9) 
given by V e Y 1 (9)(V e X 1 (9))~ 1 o-(.) (see [IS]). in the case where b, f and g are Cl w.r.t. x, 
y and z, we obtain from (I3.10p 

\Zl(6)\ < e i2La+K f )T (K g + TK f )K a a.s. 

since \(V e X 1 (9))~ 1 \ < e LaT (the proof of this inequality is similar to the one of (I3.4p ). 
When b, f and g are not differentiate, we can also prove the result by a standard approxi- 
mation and stability results for BSDEs with linear growth. □ 

Proposition 3.3. Suppose that (HF), (HBQ) and (HBQD) hold. Then, there exists a 
version of Z° such that 

| Z 0| < e 2{K f +La)T (K g + K f T) (1 + TK f e K f T (l + L a e LaT ))K a , < t < T . 

Proof. Thanks to the Malliavin calculus, it is classical to show that a version of Z° is given 
by VF (VX°)"V(.) (see [IS]). So, in the case where b, (3, g and / are Cl w.r.t. x, y, z and 
u, we obtain from (I3.3P and Lemma [3.51 

| Z 0| < e 2{K f +La)T (K g + K f T)(l + TK f e K f T (l + L a e LaT ))M a a.s. 

since KVXj )" 1 ! < e LaT (the proof of this inequality is similar to the one of (13. 3p ). 

When b, (3, g and / are not differentiate, we can also prove the result by a standard 
approximation and stability results for BSDEs with linear growth. □ 
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4 Discrete-time approximation for the FBSDE 



4.1 Discrete-time scheme for the FBSDE 

Throughout the sequel, we consider a discretization grid n := {t , . . . ,t n } of [0,T] with 
= to < ti < . . . < t n — T. For t G [0, T], we denote by ix{t) the largest element of n 
smaller than t 

7r(£) := max { ti , i = 0, . . . , n \ U < t } . 

We also denote by 1 7r | the mesh of tt 

\tt\ := max { t i+ i — ti , i = 0, . . . , n — 1 } , 

that we suppose satisfying \ir\ < 1, and by AW* (resp. At^") the increment of W (resp. the 
difference) between ^ and tj_i: AW/ 1 " := W tj - (resp. AtJ := U - for 1 < z < n. 

We introduce an approximation of the process X based on the discretization of the 
processes X° and X 1 . 

• Euler scheme for X°. We consider the classical scheme X 0,7T defined by 
X°' n - x 

~ C4 15) 

XST = + KU-i, XgjAtf + aiti^AW: , l<i<n. 

• Euler scheme for X . Since the process X 1 depends on two parameters t and 8, we 
introduce a discretization of X 1 in these two variables. We then consider the following 
scheme^ 

X%*(t k ) = x + f3(t , x)5 k=0 , < k < N , 

X^(t k ) = (t k ) + b(ti_ u X^ (t k ))Ati + aiU^AW? (4.16) 
+ P(U- X , X^ (t k ))5 i=k , 1 < z, k < n . 

We are now able to provide an approximation of the process X solution to the FSDE (12. ip . 
We consider the scheme X n defined by 

XI = X^ t) l t<T + X^{T))t t > T , 0<t<T. (4.17) 

We shall denote by {J~f ,7r }o<i< n (resp. {^ 1,7r (^)}o<i<n) the discrete-time filtration associated 
with X°> w (resp. X^{6)) 

^ ■= o-(X?!\j<i) 
(resp. J±«(0) := a(X^(9), j < i)) . 

We introduce an approximation of (Y, Z) based on the discretization of (Y°, Z°) and 
(Y 1 ^ 1 ). To this end we introduce the backward implicit schemes on 7r associated with 



1 S{ = k — if i k and 5i=k = 1 if z = fc. 
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the BSDEs (j2.4p and (12. 6p . Since the system is recursively coupled, we first introduce the 
scheme associated with (12 .4p . We then use it to define the scheme associated with (12.61 ). 

• Backward Euler scheme for (Y , Z 1 ). We consider the classical implicit scheme (Y 1,7T , Z 1 ' 71 ) 
defined by 

' Y^(n(6)) = g(X^(n(9))) , 

= E [tW«))|^i(^))] 

• Backward Euler scheme for (Y°,Z°). Since the generator of (I2.6P involves the process 
(Y^(t)) te T T j, we consider a discretization based on Y 1,7T . We therefore consider the scheme 
(yo,w z o,^ defined by 

< = eK-i^ + t^x-,^,^-)^, (419) 

where x, y, z) := f(t, x, y, z, Y^(t)) - y) for all (t, x, Me[0,T]xlxlx K. 
We then consider the following scheme defined by 

XT = ^)W + ^ t )(^(r))l f > T , 

Z* = Z°^ t) t t < T + Z^(T))l t>r , (4.20) 
Ut = (Y^(7r(t))-Y w %)t t < T , 

for t E [0,T]. 

4.2 Convergence of the scheme for the FBSDE 

We now concentrate on the error approximation of the processes X, Y, Z and U by their 
scheme X n , Y n , Z n and U n . To this end we introduce extra assumptions on the regularity 
of the forward coefficients w.r.t. the time variable t. 

(HFD) There exists a constant K t such that the functions b, a and (3 satisfy 

\b(t,x)-b(t',x)\ + \a(t)-o-(t')\ < K t \t-t'\% , 
\(3(t,x)-P(t',x)\ + \a(t)-a(t')\ < K t \t-t'\, 

for all (t,t',x) E [0,T] 2 x R. 

We can now state the main result of the paper. 
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Theorem 4.1. Under (HF), (HFD), (HBQ) and (HBQD) we have the following esti- 
mate 



E 


sup \X t -X?\ 2 


+ sup E 


\Y t -Y« 




L te[o,T] 


te[o,T] 







r f T 




r f T 




+E 


/ \z t -z?\ 2 dt 


+ E 


/ \ t \U t -U?\ 2 dt 


< K\tt\ 













for a constant K which does not depend on it. 
Proof. Fix M6R such that 

M > max {e {2La+K f )T (K g + TK f )K a ; 

e 2(K f+ L a )T {Kg + KfT) ^ + TKfe K fT{} + ^L^) ^ j ^ 

and define the function / by 

f(t,x,y,z,u) = f(t,x,y,<f M (z),u) , (t, x, y, z, u) e [0, T] x R x E x R x 

where 



* if |z| < M 
In if z > M 



We notice that <pm is Lipschitz continuous and bounded. Therefore we obtain from (HBQD) 
that / is Lipschitz continuous. 

Moreover, using Propositions I3.2l and l3.3[ we get that under (HF), (HBQ) and (HBQD), 
(X, Y, Z) is also solution to the Lipschitz FBSDE 



X f 



x 



b(s,X s )ds+ / a(s)dW s + / (3(s,X s -)dH s , 0<t<T 



Y t = g(X T )+ f(s,X s ,Y sl Z s ,U s {l-H s ))ds 



T 



Z s dW s 



U K dH„ , < t < T . 



Applying Theorem 5.1 in [16], we get the result 



□ 
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